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ON THE ORIGIN OF THE TYPE lA SUPERNOVA WIDTH-LUMINOSITY RELATION 
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ABSTRACT 

Brighter Type la supernovae (SNe la) have broader, more slowly declining _B-band light curves than 
dimmer SNe la. We study the physical origin of this width-luminosity relation (WLR) using detailed 
radiative transfer calculations of Chandrasekhar mass SN la models. We find that the luminosity 
dependence of the diffusion time (emphasized in previous studies) is in fact of secondary relevance in 
understanding the model WLR. Instead, the essential physics involves the luminosity dependence of 
the spectroscopic/ color evolution of SNe la. Following maximum- light, the SN colors are increasingly 
affected by the development of numerous Fe II/Co II lines which blanket the B-band and, at the same 
time, increase the emissivity at longer wavelengths. Because dimmer SNe la are generally cooler, they 
experience an earlier onset of Fe HI to Fe II recombination in the iron-group rich layers of ejecta, 
resulting in a more rapid evolution of the SN colors to the red. The faster _B-band decline rate of 
dimmer SNe la thus reflects their faster ionization evolution. 
Subject headings: radiative transfer - supernovae; cosmology 



1. INTRODUCTION 

Although normal Type la supernovae (SNe la) are 
generally considered a homogeneous class, they never- 
theless show substantial (~ 1 mag) variations in lumi- 
nosity at peak. The use of SNe la for cosmology mea- 
surements thus relies on empirical calibration techniques. 
Most common among the se is the width-luminosity rela- 
tion ( Phill ips et al.lll99^ . Bright SNe la generally have 
broad B-band light curves (LCs) which decline slowly af- 
ter peak. Dimmer SNe la have narrower, more quickly 
declining LCs. 

Understanding the physical origin of the width- 
luminosity relation (WLR) is a primary goal of the theory 
and modeling of SNe la. Unfortunately, a full theoreti- 
cal description of the SN la explosion mechanism is still 
lacking. Normal SNe la are widely believed to be the 
thermonuclear disruption of carbon-oxygen white dwarfs 
near the Chandrasekhar limit, however a number of un- 
certainties remain regarding the structure of the progen- 
itor, the precise ignition conditions, and the physics of 
the turbulent nuclear combustion that unbinds the star. 

Fortunately, these uncertainties need not be fully re- 
solved in order to study the WLR. In SN la explosions, 
hydrodynamic and nuclear burning processes last only 

1 minute. The subsequent luminosity is powered en- 
tirely by the decay of radioactive elements synthesized 
in the explosion, in particular ^^Ni in the decay chain 
56Ni ^ 56^0 ^ 56pg_ rpj^g j^^g gj^g therefore 

fully determined by the composition and density struc- 
ture of the material ejected in the explosion. Naturally, 
the mass of ^^Ni produced (Mni) is the primary deter- 
minate of the peak brightness of the event. Observations 
indicate that for normal objects, Mni spans the range 
0.4 - 0.9 Mq, with a typical value Mm ~ 0.6 M©. 

The challenge thus falls to radiative transfer the- 
ory to explain why SNe la with larger Mni also have 
broader B-band LCs. Most previous transfer studies 
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have e mphasized diffusion arguments of one sort or an- 
other lIHoflich et abllTgM 12001 IPinto fc Ea,stma,nll2nnib 
iMazzali et al.ll2(¥)]l) . In SNe la, the ejecta remain op- 
tically thick for the first several months after explosion. 
The width of the bolometric LC is related to the time 
scale for photons to escape the ejecta by diffusion. In 
principle, the WLR could be explained if brighter SNe la 
have higher effective opacities and hence a longer diffu- 
sion time. Different authors have invoked different physi- 
cal arguments to motivate this sort of opacity dependence 
(see 

In this paper, we stress that the WLR does not in fact 
hinge on the bolometric diffusion time, but is instead 
principally a broadband phenomenon. In particular, the 
i?-band LC depends sensitively on the rate at which the 
SN colors evolve progressively redward following maxi- 
mum light. Dimmer SNe la exhibit a faster color evo- 
lution than brighter SNe la; this turns out to be a pri- 
mary reason for their relatively faster i3-band decline 
rates. A relevant theoretical explanation of the WLR 
should therefore focus on the effect Mni has on the spec- 
troscopic/ color evolution of SNe la, rather than on the 
overall diffusion time-scale. 

In what follows, we use realistic time-dependent multi- 
group radiative transfer calculations to explain the phys- 
ical origin of the WLR in a set of Chandrasekhar-mass 
SN la models. We demonstrate that the faster color 
evolution of dimmer models reflects their faster ioniza- 
tion evolution. Following maximum- light, the SN colors 
are increasingly affected by the development of numerous 
Fe II/Co II lines which blanket the i?-band and, at the 
same time, increase the emissivity at longer wavelengths. 
Because dimmer SNe la are generally cooler, they expe- 
rience an earlier onset of Fe HI to Fe II recombination in 
the iron-group rich layers of ejecta, resulting in a more 
rapid evolution of the SN colors to the red. This explains 
their faster _B-band decline rate. 

2. PREVIOUS THEORETICAL STUDIES 

Previous radiative transfer studies of SNe la have oc- 
casionally reproduced the basic trend of the WLR among 
sets of models with varying Mni. In these studies, the 
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model relation is typically explained by claiming a strong 
luminosity dependence of the diffusion time. To the ex- 
tent that the transfer calculations include the relevant 
multi-wavelength processes, color evolution effects are 
also captured in the simulations as a matter of course. 
Sometimes their relevance t o the WLR is explic itly rec- 
ognized and quantified fe.g.. lMazzali et aDl200H) . 

Using simple dimensional arguments, one can show 
that the photon diffusion time (and hence bolometric 
LC width) scales with mean opacity as td oc k}^"^ . In 
SNe la, the dominant opacities are due to electron scat- 
tering and bound-bound line transitions, the latter of 
which is enhanced by Doppler broadening in the differen- 
tially expanding ejecta. Iron group species (iron, cobalt, 
and nickel) have the largest line opacity, due to their 
complex atomic structure. The iron group line opacities 
increases sharply to the blue, due to the larger number 
of lines at shorter wavelengths. 

At least three different (though related) arguments 
have been given as to why k (and hence the diffusion 
time) should increase with Mni, each emphasizing a dif- 
fe rent physical mechan ism. 

iHoflich et alJ ll2002l see also iKhokhlov eraP l)1993j) : 
IHoflich et all Wm^ ^ emphasize the temperature depen- 
dence of the diffusion time. Because SNe la with larger 
Mni have higher temperatures, the bulk of the radiation 
energy is concentrated at shorter wavelengths. Since the 
line opacity increases sharply to the blue, the effective 
mean opacity (and hence the diffusion time) increases 
w ith Mni- 

iPinto fc EastmanI l|2001l see also iPinto fc EastmanI 
((2000a b)) emphasize the effect of the ionization state 
on the diffusion time. They argue that is largely de- 
termined by the rate at which blue/ultraviolet photons 
fluoresce to longer wavelengths via interaction with iron 
group elements, thereby enhancing photon escape. They 
suggest that fluorescence is more efficient in singly, as 
opposed to doubly ionized species. Because SNe la with 
lower are cooler and less ionized, the diffusion time 
is shorter. 

i Mazzali et al.l l|2001j) emphasize the effect of the ejecta 
composition on the diffusion time. SNe la with large 
Mni have a greater abundance of iron group elements, 
and thus presumably larger opacities. 

In fact, all three of the above mentioned diffusion ef- 
fects must be operative at some level in the LCs of SNe la. 
The magnitude of the effects (in sum) can be quantified 
by direct examination of the model bolometric LCs. If 
diffusion time arguments are indeed of prime relevance 
in explaining the WLR, it follows that the bolometric 
LCs of brighter models must be significantly broader 
than those of dimmer models. However, in the previous 
studies this bolometric behavior is not explicitly demon- 
strated using detailed multi-group transfer calculations 
and realistic opacities. 

In the following, we present a set of model calculations 
in which the shape of the bolometric LCs in fact shows 
very little direct dependence on Mni. Variations in the 
diffusion time are therefore of secondary relevance in ad- 
dressing the B-band WLR. Instead, the model WLR pri- 
marily reflects the faster spectroscopic/color evolution of 
SNe la with lower M-^\. It is therefore of fundamental im- 
portance (and a primary goal of this paper) to explain in 



detail the physics controlling the rate of color evolution 
in SNe la. 

Because the flux mean opacity itself depends on the 
spectral energy distribution of the radiation field, the 
physics controlling the color evolution of SNe la is of 
course deeply intertwined with that determining the dif- 
fusion time. In the end, however, the relative importance 
of each can be quantified by comparing the strength of 
the WLR in the broadband versus the bolometric LCs of 
SNe la. We suggest below that such an analysis may be 
very useful in furthering our understanding of SNe la. 

3. MODELS 

In a companion paper ijWooslev et al.ll2n?)^ hereafter 
Paper I), we perform an extensive parameter study of 
SN la LCs using 1-dimensional (ID) hydrodynamical 
models of the explosion. Here we focus on a restricted 
set of those models in order to demonstrate the basic 
transfer physics underlying the WLR. Each model is a 
derivative of the spherical Chandrasekhar mass explosion 
model, M070103, which has a four zone stratified com- 
positional structure. From the center out, these zones 
consist of: O.I M© of stable iron group elements; 0.7 M© 
of s^M; 0.3 Mq of intermediate mass elements (IMEs; 
silicon, sulfur, calcium, argon); and 0.3 M© of unburned 
carbon/oxygen. Mild mixing is applied at the zone in- 
terfaces. The synthetic LCs and spectral time-series of 
model M070I03 are in good agreement with observed 
normal SNe la (Paper I). 

A simple prescription is used to vary the ^^Ni abun- 
dance in model M070I03, allowing us to isolate the effect 
of Mni on the LCs. New models are constructed by re- 
ducing the abundance of ^^Ni everywhere by a fractional 
amount (viz., 10,20,30,40, and 50%) while increasing 
the abundance of IMEs by a corresponding amount. This 
amounts to a variation of Mni from 0.35 to 0.7 Mq while 
maintaining a fixed velocity distribution of ^^Ni. Because 
the energy released in burning to IMEs is comparable to 
burning to ^^Ni, the dynamics of the explosion will not 
be greatly affected by the compositional change. The 
ejecta structures bear resemblance to som e ID delayed- 
deton ation or pulsed-detonation models l|Hoflich et alJ 

Hop) . 

We compute synthetic LCs of the models using the 
multi-dimensional time-dependent Monte Carlo radiative 
transfer code SEDDNA (Kasen et al. 2006). Given a ho- 
mologously expanding SN ejecta structure, SEDQNA calcu- 
lates the emergent spectral time series at high wavelength 
resolution. Broadband LCs are then constructed by con- 
volving the spectra with the appropriate filter transmis- 
sion functions. The code includes a detailed gamma-ray 
transfer procedure to determine the rate of radioactive 
energy deposition, and a radiative equilibrium solution 
of the temperature structure. No ad-hoc inner bound- 
ary condition is employed. Ionization and excitation are 
computed assuming local thermodynamic equilibrium, 
and bound-bound line transitions are treated using the 
expansion opacity formalism and an approximate two- 
level atom approach to wavelength redistribution, assum- 
ing a constant redistribution probability e = 0.8 for all 
lines. We have used SEDQNA to recalculate the broad- 
ba nd LCs of the scal ed models used to study the WLR 
in IPinto fc EastmanI l|200Ll and found exc ellent agree- 
ment with the results presented there. See iKasen et alJ 
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TABLE 1 

Observable properties of the models 



MNi 


tbol 


tB 


A4oi 


Mb 


Mv 


AMi5(bol) 


AA/i5(B) 




AMi5((7By/?/) 


{B — V)max 


{B - V)i5 


0.70 


15.56 


19.34 


-19.46 


-19.34 


-19.35 


0.81 


1.09 


0.72 


0.83 


-0.017 


0.568 


0.63 


15.55 


19.00 


-19.36 


-19.25 


-19.27 


0.82 


1.17 


0.77 


0.84 


-0.007 


0.621 


0.56 


15.54 


18.70 


-19.24 


-19.14 


-19.18 


0.83 


1.28 


0.82 


0.88 


0.011 


0.693 


0.49 


15.56 


18.36 


-19.11 


-19.02 


-19.08 


0.84 


1.38 


0.86 


0.91 


0.034 


0.750 


0.42 


15.52 


17.85 


-18.96 


-18.87 


-18.97 


0.84 


1.45 


0.92 


0.94 


0.057 


0.802 


0.35 


15.51 


17.22 


-18.78 


-18.71 


-18.83 


0.84 


1.53 


0.98 


0.96 


0.089 


0.851 



-19.5 




Days Since Explosion A M,^ (B) 

Fig. 1. — _B-band light curves (left panel) and width- 
luminosity relation (right panel) of the models discussed in the 
text. The shaded region in the right panel is the empirical WLR of 
[PhiUiDS ct al. 1 1999) with a magnitude dispersion of 0.15 mag and 
absolute calibration Mb = —19.3 mag at AA/i5(_B) = 1.1 mag. 



1)20061) for a detailed code description and verification. 

4. ORIGIN OF THE WIDTH-LUMINOSITY RELATION 

The WLR is often quantified as a relation between 
peak _B-band magnitude, M^, and the drop in _B-band 
magnitude 15 days after peak, AMi5(i?). Table 1 lists 
the relevant observable parameters of the synthetic LC 
models discussed in the previous section. In the mod- 
els, the variation of A/ni from 0.35 to 0.70 M© results 
in a variation in Mb from -18.71 to -19.34 mag, and in 
AMi5(B) from 1.53 to 1.09 mag. As seen in Figure 1, 
the models obey a clear WLR, with a slope that agrees 
reasonably (though not exactly) with the observed one 
of Phillips ct al. (1999). 

The model WLR obtained in Figure Q is not at- 
tributable to a strong systematic dependence of the diffu- 
sion time on Mni, as emphasized in previous studies. On 
the contrary, the bolometric rise time to maximum and 
the bolometric decline rate are similar among all models: 
tboi = 15.51 - 15.56 days, AMi5(bol) = 0.81 - 0.84 mag. 
As seen in FigureEl differences in the bolometric LCs are 
maximal (though still modest) when measured 30 days 
after B-maximum: AM3o(bol) = 1.43 — 1.61 mag. Thus, 
in this set of models, the luminosity dependence of the 
diffusion time is in fact of secondary relevance in under- 
standing the WLR. 

Instead, we identify the crucial physical effect to be 



that of Mni on the spectroscopic/color evolution of the 
models. As seen in Figure|3 all models have very similar 
colors at i?-band maximum, with B~V and B — R within 
and 0.15 mag. Thereafter, the color curves diverge, 
with the dimmer models evolving redward more rapidly. 
By 15 days after _B-maximum, differences m B — V and 
B — R color have grown to nearly 0.3 and 0.6 mag re- 
spectively. It is this faster color evolution of the dimmer 
models that accounts for their larger AMi^{B). 

The origin of this color evolution behavior can be 
traced directly to the impact of iron-group lines on the 
post- maximum spectra. Soon after maximum light, the 
i?-band is increasingly affected by the development of 
Fe II/Co II lines (Figure0J. Especially prominent are the 
strong absorption blends at wavelengths 4000-4600 A and 
4800-5300 A, however a pseudo-continuum of numerous 
weaker lines also contribute to the overall line blanket- 
ing in the blue. At the same time, a large number of 
blended Fe II/Co II lines emerge in the far red and near- 
infrared, enhancing the pseudo-continuum emissivity at 
these wavelengths (Kascn 2006). The net effect of the 
onset of the Fe II/Co II lines is thus to redistribute flux 
absorbed in the blue part of the spectrum to longer wave- 
lengths. Naturally, the rate at which these lines develop 
plays a critical role in determining the color evolution 
(and hence i?-band decline rate) of SNe la. 

In general, Fe II/Co II lines become prominent once the 
layers of ejecta rich in ^^Ni (since decayed to iron/cobalt) 
begin to recombine from doubly to singly ionized, a pro- 
cess which occurs quite suddenly at temperature T « 
7000 K. For the epochs at and before i?-band maximum, 
the iron/cobalt-rich layers of ejecta [v < 9000 km s""'^) 
have T > 7000 K and Fe II/Co II lines are weak or ab- 
sent in all model spectra (see Figure |S1 top panel). At 
these times, the optical continuum is pseudo-blackbody 
and the spectral and color differences among the models 
are minor, even when A/ni is varied by a factor of two. 
The brighter/broader SNe la have a bluer {B — T^)max 
color, but by a relatively smal l amount, as noted in ob- 
servations (IPhillins et a,].llT99 l. 

In the weeks following maximum light, however, the 
layers of iron-rich ejecta progressively cool to T < 7000 K 
and, as these layers recombine, the Fe II/Co II lines be- 
come steadily stronger. The exact rate at which the lines 
develop depends sensitively on the overall ejecta temper- 
ature scale. Models with lower Mni have generally lower 
temperatures and hence an earlier onset of Fe III to Fe II 
recombination. As seen in Figure [S] (bottom panel) by 15 
days after maximum the line blanketing is much stronger 
(and the colors much redder) in the lower A^ni model. 
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The larger AA'/i5(i?) of the dimmer model thus reflects 
the more rapid development of the Fe II/Co II lines. 

To be clear, the WLR is not merely a consequence of 
the generally redder colors of dimmer SNe la, but rather 
of the faster rate at which the colors evolve in these ob- 
jects. In considering this spectroscopic effect, it is useful 
to keep in mind that the temperature differences among 
normal SNe la are expected to be rather small. The 
ejecta temperature (as determined by the balance of ra- 
diative heating and cooling) scales with luminosity as 

T oc L^/^ and hence T cx M^(^. A factor of two differ- 
ence in Mni therefore amounts to only a ~ 20% differ- 
ence in temperature. If the SN spectrum is reasonably 
characterized by a blackbody, the corresponding varia- 
tions in B — V or B — R color would be relatively small 
(< 0.15 mag). Indeed, this is essentially why the pre- 
maximum and maximum light colors of normal SNe la 
show only a mild dependence on Mni- 

In the post-maximum epochs, however, the SN la 
spectrum increasingly deviates from blackbody due to 
the burgeoning strength of the Fe II/Co II lines. The 
strength of these lines depends very sensitively on the 
ionization state. The ionization in turn obeys an ex- 
tremely non-linear dependence on temperature, under- 
going sudden recombination at T « 7000 K. Thus, in the 
post-maximum epochs, small differences in temperature 
translate to large differences in B — V and B — R color. 
The onset of this non-linear behavior is why the decline 
rate /S.Miz{B) offers a good diagnostic of the SN tem- 
perature scale, and hence luminosity. Other approaches 
to calibratiiig SN e la, such as the CMA GIC method 
(iWang et al.ll2003|) and the C12 correlation ijWang et all 
12005ft are likewise explained by the same strong depen- 
dence of the SN colors on the ionization state during the 
weeks following maximum light. 

At around 30 days after _B-maximum, the iron group 
layers of ejecta finish completely recombining from dou- 
bly to singly ionized. From this time on, differences in 
the ionization state among the models cease to be of great 
significance. The color curves reconverge and the mod- 
els enter a phase of nearly identical color evolution (see 
Figure ^ . The s ame g eneral behavior is noted in the 
observations ijLiral 119961 IPhillips et al.lll999j) . For this 
reason, the decline rate measured over longer time peri- 
ods (e.g., AM3o(i?)) provides a less sensitive measure of 
the model peak magnitude. 

5. RELEVANCE TO OBSERVATIONS 

The observed color evolution of SNe la bears a strik- 
ing similarity to the model behavior of Figure (e.g., 
I^ang ct al. 2005, Figure 1). Normal SNe la display fairly 
uniform B — V colors at maximum light; thereafter the 
dimmer objects become redder more quickly. The ob- 
served differences in color at day -1-15 are comparable to 
the differences in AMi5(i?). 

FolatcUi (200^ quantifies the spectral evolution of 
SNe la directly using equivalent width measures. His ap- 
proach nicely illustrates the sudden growth of Fe II/Co II 
features (there labeled "Mg 11") in the B-band begin- 
ning just after maximum light. The rate at which 
these features develop is shown to correlate tightly with 
AMi5(i?). Such analyses clearly point to the spectro- 
scopic nature of the WLR. 

The construction of bolometric LCs from uvoir obser- 





-20 20 40 60 80 
Days Since Bolometric Maximum 



-20 -10 10 20 30 40 50 
Days Since B-Band Maximum 



Fig. 2. — Left: Bolometric LCs of the models, normalized at peak 
and shown relative to time of bolometric maximum. The LC shapes 
are all very similar over the first three weeks after maximum light. 
Moderate differences of order 0.15 mag appear around 30 days 
after maximum, when the brighter models reach a small secondary 
maximum. Right: B-band of the models, normalized at peak and 
shown relative to time of B-maximum. The LC decline rate just 
after maximum depends strongly on Mni. 
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Fig. 3. — Color evolution of the models discussed in the text. In 
all panels, the dimmer models are the topmost lines. All models 
have similar colors near maximum light; thereafter the dimmer 
ones (lower M^i) become redder more rapidly. The differing rates 
of color evolution account for differences in AAfi5(B). 



vations of SNe la can be used to further explore this issue. 
If the luminosity dependence of the diffusion time is in 
fact central to explaining the WLR, one expects to ob- 
serve a strong width luminosity relation in the bolometric 
LCs of SNe la. It is currently unclear from the limited 
observational sample whether such a relation exists at all. 
Phillios ct al. (2006), for example, construct bolometric 
LCs from U BVRI- and (when available) JHK-ha,nd ob- 
servations of 15 SNe la. The bolometric decline rates 
AMi5(bol) of the sample show some degree of scatter, 
but no clear correlation with the peak bolometric mag- 
nitude Afboi- 

The bolometric LC constructions of ' Contardo et alJ 
(^000), which use primarily i?V^i?/-band observations, 
do show a noisy correlation between AMi5(bol) and 
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Fig. 4. — Spectral evolution from S-band maximum to twenty 
days later for the Mni = 0.7 Mq model. The figure illustrates how 
the drop in _B-band magnitude depends sensitively on color changes 
due, in large part, to the development of spectral line blanketing 
from iron group elements. The shaded bands in the plot highlight 
two prominent spectral absorption features due to Fe 11/ Co II line 
blends. 
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Fig. 5. — Demonstration of the crucial effect of Fe II/Co II 
line blanketing on AMir^{B). The figure compares synthetic spec- 
tra of the AfNi= 0.35 and 4/^1= 0.70 Mq models. All spectra 
have been normalized by dividing by the bolometric luminosity at 
the time. The blue, green, and red color shadings show the re- 
gions of the B,V and R bands respectively. Near maximum light 
(top panel) both models have similar, pseudo-blackbody coninuua. 
Fifteen days later (bottom panel) broad absorption features from 
Fe II/Co II lines dominate the B-band. Because the cooler MNi= 
0.35 M0 model is less ionized, the post-maximum line blanketing 
develops more quickly, leading to the faster B-band decline rate. 



Mbol (see also IStritzinger et al.ll2006|) . This correlation, 
however, may in large part reflect the lack of ultravio- 
let and near-infrared (NIR) data used in the bolometric 
constructions. At day 15 after i?-maximum, for exam- 
ple, 20% or m ore of the SN la flux may occur at NIR 
wavelengths (Suntzcfl 2003). As the NIR contribution 
at day 15 is relatively greater in the dimmer SNe la, 
UBVRI constructions may systematically overestimate 



the value of AAfi5(bol) for dimmer objects by as much 
as 0.1 — 0.2 mag. This behavior is seen in our models 
(Table 1) where one notes that AMi5{UBVRI) exhibits 
a much stronger dependence on A/ni than does the true 
bolometric decline rate AA/i5(bol). In fact, the slope 
of our UBVRI model WLR agrees fai rly well with the 
empirical one of iContardo et alJ l)200(H) . suggesting that 
the observed correlation reflects, in large part, the role 
of color evolution. 

Clearly the construction of true bolometric LCs for a 
sample of SNe la using UV through X-band observa- 
tions would be of great value in constraining the radiative 
transfer in SNe la. In the speciflc models studied in this 
paper, differences in the bolometric LCs are relatively 
small, and attributable to the reduction in flux mean 
opacit y which accompan i es the onset o f Fe II/Co II lines 
(Pin to fc Eastmanll200lt lKasenll200d|) . The bolometric 
differences become most prominent around 30 days af- 
ter i3-maximum, when the brighter models show a slight 
rise to a secondary bolometric maximum. Differences in 
diffusion time may therefore be somewhat more relevant 
in formulations of the WLR employing later time data. 

The bolometric LCs of SNe la may also depend sensi- 
tively on the ejecta kinetic energy and the radial distri- 
bution of ^^Ni (see Paper I). By varying these two other 
parameters, it is possible to construct sets of models that 
show greater or lesser degrees of scatter and/or correla- 
tion in the bolometric WLR. In addition, NLTE effects 
(not included here) may also contribute to a stronger lu- 
minosity dependence of the bolometric LCs, especially 
if the probability of wavelength redistribution in lines is 
strongly temperature dependent. Such issues will be the 
subject of future studies. However, because the rate of 
color evolution dominates the B-band decline rate, the 
WLR should remain a robust luminosity calibration even 
when the bolometric LCs of SNe la show significant ran- 
dom variations. 

The bolometric rise time to maximum is the same for 
all models in this paper: tboi = 15.5 days. The S-band 
rise time (^3)7 on the other hand, depends on the color 
evolution. Because iron group recombination occurs ear- 
lier in lower Mm models, their colors begin evolving red- 
ward at earlier epochs. Thus, in both models and obser- 
vations, the dimmer SNe la exhibit a shorter te- 

Because Fe II/Co II line blanketing generally increases 
to the blue, the WLR will usually be stronger in bluer 
wavelength bands. The V-band, for instance, is impacted 
by the development of Fe II/Co II lines to a lesser extent 
than the B-band. This leads, in both models and obser- 
vations, to a slower F-band decline rate AMi5{V) and 
a weaker dependence of AAfi5(y) on the peak T^-band 
magnitude Afy. 

The i?-band decline rate also has an intimate connec- 
tion to the double-peaked morphology of the far red and 
NIR LCs of SNe la. The flux absorbed in Fe II/Co II lines 
at blue wavelengths is redistributed to the red, eventually 
leading to a NIR secondary maximum ( Pi nto fc Eastmanl 
.2Q00b; Kasen 2QD,6) . Naturally, in both models and ob- 
servations, SNe la with faster B-band decline rates have 
earlier NIR secondary maxima. 

On the whole, one recognizes that the ionization evolu- 
tion of iron group elements is the most significant physi- 
cal factor in understanding SN la LCs and the WLR. The 
faster overall spectroscopic and photometric evolution of 
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dimmer SNe la reflects, primarily, their faster ionization 
evolution. 

In addition, the distribution of iron group elements in 
the ejecta can be identified as a crucial secondary param- 
eter in the WLR. All models in this paper have substan- 
tial iron group abundance out to v 9000 km s~^. If 
the iron group elements are mixed to yet higher veloci- 
ties, the development of Fe II/Co II hues will be hastened 
and intensified. Thus, for given Mnj, enhancing the iron 
group abundance in the outer layers of ejecta increases 
the i?-band decline rate. The relatively low dispersion in 
the observed WLR suggests that normal SNe la exhibit 
limited variability in their iron group distributions. The 
physical origin of this regularity, however, can only be 
addressed by detailed first principle explosion modeling 
of SNe la. 

Fortunately for future precision cosmology experi- 
ments, the radiative transfer effects identified in this 
paper can be used to monitor and limit deviations in 
the WLR. Variations in the velocity distribution of iron 
group elements, for example, can be constrained by mea- 
suring the Doppler-shifts of iron features in the SN la 
spectra. NIR observations provide additional indicators 
of the ioniz ation evolutio n and iron group distribution in 
the ejecta (IKasenl|2 006!). Close inspection of well sam- 
pled color curves (or spectral time series) should reveal 
equivalent or improved calibration techniques and may 



be useful in quantifying subtle deviations from the stan- 
dard evolution. These and related tests should be feasible 
with the data from ongoing observational campaigns, and 
will help to limit the intrinsic scatter or evolutionary bi- 
ases potentially affecting SNe la cosmology experiments. 
With the aid of detailed transfer models such as those 
presented here, refinements to the cahbration techniques 
need not rely on the fortuitous discovery of empirical 
correlations, but can be devised to reflect the essential 
physics of SN la light curves. 
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